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A  theoretical  investigation  was  made  of  the  non-equilibrium  expansion  of  a  par¬ 
tially  ionized  gas  through  a  supersonic  nozzle.  Both  hydrogen  and  argon  were  studied. 
The  ionization  and  recombination  rate  parameters  used  were  those  calculated  by  Bates 
and  co-workers,  using  their  collisional-radiative  model  of  the  recombination  process. 
These  calculations,  for  hydrogen,  include  the  influence  of  radiation  trapping  on  the 
overall  rates,  and  yield  as  well  the  amount  of  recombination  energy  which  is  gained 
by  the  third-body  electron  in  the  recombination. 

The  energy  balance  for  the  electrons  and  massive  particles  was  studied  for  both 
the  optically  thin  (all  recombination  radiation  lost)  and  optically  thick  (all  radia¬ 
tion  absorbed)  cases.  It  was  found  that  the  recombination  process  produces  an 
increase  in  electron  temperature  over  that  of  the  ion-atom  temperature,  and  that 
this  temperature  difference  is  greater  for  the  optically  thick  than  for  the  optically 
thin  case.  Also,  there  is  more  net  recombination  for  the  optically  thin  case. 

For  high  initial  ionization  (62%),  recombination  affects  the  gross  flow  variables 
measurably,  when  compared  with  the  "frozen  flow"  solution.  For  low  initial  ionization 
(1%)  the  gross  flow  variables  are  practically  unaltered,  and  the  principal  effect  of 
the  recombination  is  to  produce  an  elevated  electron  temperature. 

On  comparing  the  hydrogen  and  argon  flows,  one  finds  that  the  effect  of  the 
larger  atomic  mass  of  the  argon  is  to  diminish  the  effectiveness  of  electron-atom 
energy  exchange,  resulting  in  a  larger  difference  between  electron  and  atom  temper¬ 
atures.  For  large  initial  ionization,  the  fractional  amount  of  recombination  in  the 
supersonic  part  of  the  nozzle  was  greater  in  argon  than  in  hydrogen,  due  to  the  longer 
residence  time  of  the  argon  iD  the  nozzle.  But  for  low  initial  ionization,  the 
reverse  was  true,  because  the  higher  electron  temperature  in  the  argon  reduced  the 
recombination  rate. 

Electron  thermal  conductivity  effects  were  found  to  be  unimportant  for  the 
frozen  and  non-equilibrium  flows,  but  should  be  included  in  the  equilibrium  flow 
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calculations 
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1.0  INTRODUCTION 

The  purpose  of  the  present  paper  is  to  present  a  detailed  study  of  the 
flow  of  a  partially-ionized  gas  through  a  convergent-divergent  nozzle,  taking 
into  account  electron-ion  recombination  and  the  radiative  and  collisionally- 
transferred  energy  liberated  in  the  recombination.  Considerable  recent  work 
has  been  done  £ 1 -7 3  which  establishes  the  fact  that  in  many  recombining 
plasmas  the  recombination  mechanism  is  dominantly  electron-electron- ion 
three  body  recombination  to  a  highly  excited  atomic  level,  followed  by 
collisional  and  radiative  de-excitation  to  the  ground  state.  As  a  consequence 
of  this  process,  one  finds  that  a  part  of  the  potential  energy  of  ionization 
is  given  to  the  third  body  electron,  with  the  result  that  the  temperature  of 
the  electron  gas  is  elevated  above  that  of  neutral  gas.  The  energy  liberated 
in  the  form  of  radiation  may  or  may  not  be  re-absorbed  by  the  neutral  atoms, 
depending  on  the  character  of  the  plasma.  An  important  feature  of  the  re¬ 
combination  process  is  that  the  net  recombination  rate  is  strongly  dependent 
on  the  electron  temperature. 

The  collisional-radiative  recombination  model  has  been  applied  to  the 
case  of  a  decaying  static  plasma,  and  found  to  be  in  excellent  agreement 
with  experiment.  [l,3]  However,  the  consequences  of  the  model  have  not  yet 
been  examined  fully  for  a  flowing  plasma,  in  which  there  are  strong  coupling 
effects  between  the  fluid  dynamical  phenomena  and  the  recombination  parameters. 
The  purpose  of  the  present  paper  is  to  investigate  these  effects.  Ue  have 
chosen  to  examine  the  consequences  of  the  collisional-radiative  recombination 
process  in  the  rapid  expansion  of  a  partially  ionized  plasma  through  a  super¬ 
sonic  nozzle.  In  particular,  we  are  interested  in  how  the  inequality  between 
electron  and  ion  temperatures,  which  is  a  crucial  factor  in  the  recombination 
process,  differs  in  the  flow  case  from  what  is  observed  in  a  decaying  static 
plasma.  Also,  we  wish  to  ascertain  how  the  flow  is  affected  by  the  presence 
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or  absence  of  radiation  trapping.  We  have  chosen  hydrogen  as  our  model  gas, 
because  the  recombination  theory  is  most  complete  in  this  case.  However,  we 
have  also  carried  out  approximate  calculations  for  argon,  in  order  to 
ascertain  the  influence  of  the  atomic  mass  on  the  process. 

The  expansion  of  an  ionized  gas  through  a  supersonic  nozzle  has  been 
investigated  by  Bray  and  Wilson  [8],  Sockel  [9],  and  others,  without,  however, 
including  radiation  losses.  Bray  and  Wilson  studied  argon,  and  used  a  re¬ 
combination  rate  coefficient  deduced  from  the  ionization  rate  given  by 
Petschek  and  Byron  [lO].  This  procedure  is  most  appropriate  for  a  highly 
ionized  gas.  Sockel  used  an  approximate  recombination  rate  formula  taken 
from  Hinnov  and  Hirschberg  [l],  which  is  valid  for  low  electron  temperatures 
(kTg  $  0.25  eV).  Sockel  also  assumed  the  electron  and  atom- ion  temperatures 
to  be  equal.  Another  study  closely  related  to  the  present  work 
is  that  of  Byron,  Bortz  and  Russell  [l9].  These  authors  have  discussed  in 
detail  the  important  radiative  processes  and  their  influence  on  the  collisional- 
radiative  recombination  rate,  and  have  carried  out  several  calculations  on 
the  flow  of  helium  through  a  supersonic  nozzle.  We  shall  have  more  to  say 
of  this  work  later. 

In  the  present  work,  values  for  the  rate  coefficients  and  radiation 
power  are  taken  from  the  calculations  of  Bates,  Kingston  and  McWhirter  [4,5]. 

The  general  theory  of  these  calculations  can  be  found  in  the  original  papers, 
and  will  not  be  discussed  here.  For  simplicity,  two  extreme  cases  regarding 
radiative  power  are  considered.  The  first  is  the  "optically-thin"  case, 
denoted  by  0(A)  *  0,  where  a (A)  is  the  optical  thickness  of  the  plasma 
at  a  wavelength  X.  In  this  case,  all  the  radiated  energy  is  lost  from  the 
plasma.  The  other  extreme  is  the  optically  thick  case,  *(X)  ■  «,  in  which 
all  the  radiation  is  re-absorbed  by  the  plasma.  Generally  speaking,  cases 
of  practical  interest  such  as  the  plasma-jet  flows  fall  somewhere  between 
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these  two  extremes,  as  is  pointed  cut  by  Byron,  ec.«u.  119J.  A  typical  situa¬ 
tion  [19]  is  where  the  plasma  is  optically  thin  to  continuum  radiation  and 
transitions  between  excited  states,  but  optically  thick  for  transitions  to 
the  ground  state.  Bates,  ej:.al_. ,  have  in  fact  computed  rate  constants  for 
this  case,  as  well. 

Throughout  the  analysis,  the  following  assumptions  are  employed:  The 
flow  is  taken  to  be  quasi-one-dimensional.  The  self -relaxation  times  (the 
times  for  randomization  of  the  thermal  velocities)  of  all  the  species  are 
much  shorter  than  the  characteristic  flow  time,  so  that  the  velocity  distri¬ 
bution  for  each  specie  is  Maxwellian  (though  not  necessarily  at  the  same 
temperature).  Viscosity,  mass  diffusion  and  heat  conduction  effects  are 
assumed  to  be  absent.  Axial  heat  conduction  by  the  electrons  is  shown  to 
be  negligible  by  numerical  examination  of  the  final  results.  We  also  assume 
that  the  mass  flow  velocities  of  electrons,  ions  and  neutral  atoms  are  all 
equal,  and  that  the  ion  and  atom  temperatures  are  the  same.  The  equality 
of  mass  flow  velocities  implies  that  the  net  electrical  current  is  zero, 
and  that  as  a  consequence  an  electric  field  between  the  ions  and  electrons 
must  exist.  This  electric  field  plays  a  significant  role  in  the  electron 
energy  balance,  and  acts  to  equalize  the  temperatures  between  the  ions  and 
electrons. 

The  mean  free  paths  for  collisions  between  species  are  all  assumed  to 
be  much  smaller  than  the  characteristic  flow  length  (say,  the  nozzle  throat 
diameter)  so  that  continuum  flow  equations  are  applicable.  Lateral  diffusion 
and  wall  recombination  are  neglected,  and  it  is  assumed  that  no  applied  mag¬ 
netic  or  electric  fields  exist.  Also,  the  ideal  gas  equation  of  state  is 
taken  for  all  species,  and  the  additive  law  of  partial  pressures  is  employed. 

Equilibrium  solutions  based  on  the  Saha  equation  for  electron  density, 
and  "frozen-flow"  solutions  are  also  presented  for  comparison  with  the 
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finite-rate  calculations.  Strictly  speaking,  one  can  have  neither  equilibrium 
nor  frozen  flow  for  the  0*0  case,  because  in  this  case  the  flow  is  not 
adiabatic.  However,  in  the  case  o  *  «  the  equilibrium  and  frozen  flow 
solutions  should  represent  the  two  limiting  cases. 


2.0  GOVERNING  EQUATIONS 

We  take  u  to  be  the  flow  velocity,  x  the  distance  measured  along  the 
nozzle  axis,  whose  local  cross-sectional  area  is  A,  and  we  let  Ne,  N^,  Na 
be  the  electron,  ion  and  atom  number  densities,  respectively.  The  conserva¬ 
tion  of  ex  t  number  density  is  expressed  by 

dN  > 

A  ST  “  S  (Ne  A  u)  (2,1) 


and  the  global  continuity  equation  is 


-  (p  u  A)  -  0 


where 


p  =  N  m  +  N  m.  +  N  m  83  (N  +  N  )  m 
aa  ii  ee  a  ea 


(2.2a) 


(2.2b) 


We  have  assumed  the  ionic  and  atomic  masses  to  be  equal,  m^  *  m^,  and  have 
invoked  the  condition  of  plasma  neutrality  N^  *  Nfi.  The  global  momentum 
equation  is 

pu  &  *  -  &  (2-3a) 


where  we  take  for  the  total  pressure 


P  -  Pi  +  pa  +  pe  S  0»i  +  Ma)  W  +  Ne  kTe  ,  (2.3b) 

making  use  of  the  ideal  gas  law  approximation  with  T^  *  T&  *  T.  In  the 
above,  k  is  Boltzmann's  constant. 
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For  the  global  energy  equation  we  have 


in  which 
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R 

P 


i  - 


£ 
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(2.4a) 


(2.4b) 


is  the  enthalpy  per  unit  mass,  and  R  is  the  power  lost  by  radiation,  per 
unit  volume.  If  X  is  the  ionization  energy  per  atom,  then  the  internal 
energy  per  unit  volume  is 


-  !». 


kT  + 
e 


(N  +  N.)  kT)  + 
a  x 


N  X 
e 


and  the  enthalpy  and  pressure  can  be  written  as 


i  =  |-{T  +  fT)+-X 

2  m  v  e  m 

a  a 


(2.5a) 


p  «  jj-  k(T  +  f  T  ) 
m  e 

a 

where  the  degree  of  ionization  is  defined  as 

N 


(2.5b) 


N  +  N 
a  e 


(2.6) 


Equations  (1)  and  (2)  combine  to  give 

dN 


„  u  df 
m  dx 


dt 


(2.7) 


The  rate  equation  which  determines  dN  /dt  is,  in  the  collisional-radiative 


model  of  Bates,  Kingston  and  McWhirter  [4], 

2 

XT’  “  -  a  N  +  S  N  N 

dt  e  e  a 


(2.8) 


in  which  a  and  S  are  the  recombination  and  ionization  rate  constants.  It 
may  be  noted  that  the  term  involving  S  is  important  only  when  the  electron 
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density  is  near  Saha  equilibrium.  From  (2.7)  and  (2.8),  we  obtain 

Af  N 

2“  «  -  Ua  +  S)f  -  S]  ~  (2.9) 

One  more  equation  is  needed  to  complete  the  formulation  of  the  problem, 
an  electron  energy  balance  relationship.  This  is  given,  for  a  recombining 
plasma,  by 

»e  “  5  (‘e  +  2  “')  *  F“  '  -(f^e  +  x)(^)-«-Qel-0M 

(2.10) 

Here  F  is  the  force  exerted  on  the  ions  by  the  electrons,  due  to  the  induced 
electric  field  which  is  required  to  maintain  =  ue  =  u;  hence  Fu  is  the 
rate  at  which  “flow  work"  is  expended  by  the  electrons.  The  first  two  terms 
on  the  right  hand  side  of  (2.10)  together  represent  the  net  rate  of  energy 
gain  of  the  electron  gas  due  to  the  three-body  recombination  process,  it 
being  the  total  recombination  energy  rate  less  the  power  lost  by  radiation. 

The  quantities  and  Q^a  represent  energy  loss  rates  due  to  electron- 

ion  and  electron-atom  collisions,  respectively.  For  Qei»  we  use  the 
expression  given  in  Ref.  10, 


N  2e4  /  8x 


I  e  /  8X  m  /  T  -  T  s  r  9(.kT  )  -j 


9(kT  )' 


(2.11) 


in  which  e  is  the  electronic  charge.  If  we  assume  that  the  intermolecular 
force  between  electrons  and  atoms  follows  the  inverse  fifth-power  law,  then 


from  Ref.  11 


Q  -  2.5  —  (m  B)1/2  k  K  K  (T  -  T)  (2.12) 

ea  a  v  e  ea'e 

a 

where  B,  the  intermolecular  force  constant,  will  be  taken  from  experimental 
data.  (See  Appendix  S.) 


7. 


The  force  F  acting  on  the  electrons  can  be  evaluated  from  the  electron 
momentum  equation, 


P  U  -7- 
s  dx 


+  F  +  P,  +  P 
dx  ie  ae 


(2.13a) 


v  * 


where  and  P^e  represent  momentum  transfers  due  to  collisions  between 

electrons  and  massive  particles.  These  terms  vanish,  however,  if  the  assump¬ 
tion  u  «  u  **  u  is  made.  Now,  the  electron  pressure  will  in  general  be 
1  6  & 

of  the  same  order  as  the  gas  pressure  p,  but  Pg  u(du/dx)  will  be  of  order 

(m  /m  )p(du/dx).  Hence  to  good  approximation,  we  may  write 
c  fl 

dp 

F  -  (2.13b) 

with  neglect  of  terms  of  order  (m  /m  );  this  relation  implies  that  the 

force  acting  on  the  electrons  which  balances  the  electron  pressure  gradient 

is  the  induced  electric  field  E  set  up  between  the  ions  and  electrons  in 

order  to  maintain  u^^  *  ug.  Exactly  the  same  conclusion,  of  course,  follows 

from  a  consideration  of  Ohm's  law  (cf.  Spitzer,  Ref.  12,  p.  28),  if  the 

conduction  current  is  set  equal  to  zero  and  in  fact  F  -  Ng  e  E.  It  has 

been  pointed  out  by  Morse,  Mitchner  and  Kruger  [13]  that  in  order  for 

Eq.  (2.13b)  to  apply  the  nozzle  must  be  insulated,  otherwise  the  Induced 

electric  field  would  be  shorted  out.  Equation  (2.13b)  still  could  be  valid 

for  a  conducting  nozzle,  if  the  plasma  were  separated  from  the  nozzle  wall 

by  a  non-conducting  boundary  layer. 

On  substitution  of  Eqs.  (2.11),  (2.12)  and  (2.13b)  into  (2.10)  and 
1  2 

neglecting  j  u  as  small  compared  to  ig,  we  obtain 

dT  ,  T  ._  (T  +  |  r  )  .*  2  m  R 

_ e  2  e  dp  \  e  3  k  7  df  _  a 

dx  3p  dx  f  dx  ’  3p  f  u  k 


l  [Bj  £  T'3/2  m  (  ^)+  (1  -  f)  d3  ]  <Te  -  I) 


(2.14) 
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in  which 


v.  I 
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Vs 
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k  j 


Is, 

r«- 


2e4  /  8,1  “e 


3m  2k 
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/  as  m  \ 

(t*) 


1/2 


9k3m 


8ne 


5«(meB) 


1/2 


3tn 


(2.15) 
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The  problem  is  now  completely  formulated,  and  consists  of  the  governing 
equations  (2.2a),  (2.3a),  (2.4a),  (2.9)  and  (2.14),  together  with  the  state 
equations  (2.5).  The  quantities  a,  S  and  R  are  taken  from  the  tables 
of  Bates,  et.al.(4].  The  first  two  quantities,  for  the  case  <t  «  »,  are 
plotted  in  Fig.  1,  and  compared  with  the  values  used  by  Bray  and  Wilson  [8,15], 

Analytical  representations  for  these  quantities  are  presented  in  Appendix  A 
of  this  report. 

3.0  THE  ELECTRON  TEMPERATURE  EQUATION 

It  is  of  interest  to  form  a  differential  equation  for  the  electron-ion 
temperature  difference  (Tfi  -  T).  By  combining  Eqs.  (2.3a),  (2.4a),  (2.5) 
and  (2.14),  one  obtains 

d_/T  ,  «  .  2  ma  R  1  (  .  2X\df  2dfej  ,  ,  } 

dx  (Te  T)  3  p  u  f  k  f  V  e  3  k  )  dx  +  3  dx  Ue 


-  P 


II 


[Di£  1-3/2  fa(D25r)+<1-£)D3] 


(T  -  T) 
e 


(3.1) 


Now,  if  there  is  no  recombination,  so  that  R  and  df/dx  are  zero,  the 
formal  solution  of  (3.1)  can  be  written  down  as 


(Tfi  -  T) 


C  exp 


[  /  H 


(3.2a) 


where  C  is  a  constant  and 
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X pwv* jvvw  ? ' 


2  l^.n.Pl  .  Q  ,(L.f,  I), 


dx 


Dxf  T 


■3/2 


in 


/  D2Te  X  1 

(in)"1-**  D3j 


(3.2b) 
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From  (3.2)  we  see  that  if  Tg  *  T  initially  at  x  »  0,  then  C  a  0  and  T£ 
will  remain  equal  to  T  throughout  the  expansion.  This  conclusion  is  valid 
provided  heat  conduction  effects  are  negligible,  as  discussed  later  on.  If 
on  the  other  hand,  the  net  recombination  rate  is  not  zero,  then  the  electron 
gas  will  gain  energy  due  to  the  recombination  process,  and  Tg  will  in 
general  be  greater  than  T.  The  first  two  terms  on  the  right  of  (3.1) 
together  give  the  net  rate  of  energy  gain  of  the  electrons.  On  the  other 
hand,  the  negative  pressure  gradient  acts  to  decrease  the  temperature  differ¬ 
ence  (Tg  -  T),  as  does  collisional  energy  exchange  between  the  electrons 
and  the  massive  particles.  These  effects  are  represented  by  the  third  and 
fourth  terms  on  the  right  of  (3.2),  respectively.  Therefore,  the  resultant 
rate  of  change  of  (Tg  -  T)  represents  a  balance  between  the  energy  gain 
and  energy  loss  mechanisms,  and  is  determined  by  their  relative  effectiveness. 

The  effectiveness  of  collisional  energy  exchange  is  directly  proportional 
to  PD^  and  PD y  which  we  see  from  (2.15)  vary  inversely  with  ma<  Thus, 
the  larger  the  atomic  mass,  the  less  effective  will  be  the  energy  exchange 
between  the  electrons  and  the  ions  and  atoms,  and  we  expect  therefore  that 
the  temperature  difference  (Tg  -  T)  will  be  larger  in  an  argon  than  in  a 
hydrogen  expansion.  This  expectation  is  borne  out  by  the  calculations.  The 
pressure  gradient  term  also  acts  to  equalize  the  two  temperatures,  but  its 
effect  is  small  owing  to  the  logarithmic  dependence.  On  the  other  hand,  if 
the  expansion  is  rapid,  the  factor  p/u  multiplying  the  collisional  exchange 
term  decreases  rapidly,  diminishing  the  effectiveness  of  collisional  equali¬ 
zation  of  the  temperatures.  The  net  result  is  that  the  more  rapid  is  the 
expansion,  the  larger  is  the  temperature  difference. 


The  net  electron  energy  gain  depends  on  the  rate  of  recombination  and 
the  amount  of  energy  radiated  away.  Loss  of  recombination  energy  by  radia¬ 
tion  diminishes  the  amount  of  this  energy  which  is  given  to  the  electrons, 

thus  reducing  the  electron  temperature.  However,  the  recombination  rate 
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coefficient  increases  (0!  ~  Tg  '  ,  according  to  the  approximate  analytical 
representation  of  Hinnov  and  Hirschberg  [  1 ] ),  so  a  decrease  in  electron 
temperature  increases  the  rate  of  energy  release,  thus  tending  to  drive  the 
electron  temperature  in  the  opposite  direction.  As  we  shall  see,  radiation 
loss  actually  turns  out  to  decrease  the  temperature  difference  (Tg  -  T). 

4.0  EQUILIBRIUM  AND  FROZEN  FLOW  EQUATIONS 

Before  proceeding  to  the  finite-rate  calculations,  we  shall  discuss  the 
two  limiting  cases  of  equilibrium  and  frozen  flow.  In  equilibrium  flow,  the 
ion  and  electron  temperatures  are  everywhere  equal,  and  the  degree  of  ioniza¬ 
tion  is  given  by  the  Saha  equation,  which  is 

3/2 

i;a  -  C(mr)  [;3:]“»(s0  <4-! 


(4.1a) 


where  g  i3  the  quantum  mechanical  statistical  weight.  When  numerical  values 
for  hydrogen  are  inserted,  we  obtain 
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where  p  is  the  numerical  value  of  the  pressure  in  dyne/cm4'  and  T  is 
numerically  equal  to  the  temperature  in  °K.  Since  no  energy  is  radiated 
away,  the  energy  equation  becomes 
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Using  (4.2),  (2.5),  (2.3)  and  (4.1a),  we  arrive  at  the  rate  equation  in  the 
form 
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Equations  (4.1)  and  (4.2)  together  with  (2.2),  (2.3),  (2.5)  and  (2.6)  specify 
the  equilibrium  flow  problem. 

For  frozen  flow,  f  *  f  *  constant,  R  »  0  and  T  «*  T.  The  governing 
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equations  are  then  (2.2),  (2.3),  (2.5),  (2.6)  and  (4.2). 


5.0  METHOD  OF  SOLUTION 

If  we  completely  specify  the  nozzle  contour  in  advance,  the  frozen  and 
equilibrium  flow  solutions  can  be  calculated  directly,  but  it  is  necessary 
to  iterate  on  the  finite-rate  equations  in  order  to  find  the  value  of  the 
mass  flux  which  for  a  given  set  of  inlet  conditions  yields  sonic  flow  at 
the  nozzle  throat.  Since  the  finite-rate  equations  are  the  most  complicated, 
it  is  preferable  to  use  a  scheme  which  allows  these  equations  to  be  solved 
directly,  and  perofrm  the  iterations  on  the  equilibrium-  and  frozen-flow 
equations  instead.  This  can  be  accomplished  by  specifying,  for  the  finite- 
rate  equations,  the  pressure  distribution  in  the  convergent  portion  of  the 
nozzle,  and  computing  with  this  pressure  distribution  until  sonic  flow  is 
reached.  Beyond  this  point,  in  the  supersonic  portion  of  the  nozzle,  A(x) 
can  be  prescribed.  In  the  course  of  the  c.  .lculation,  A(x)  in  the  subsonic 
portion  of  the  nozzle  is  determined  from  the  continuity  equation.  Iterations 
on  the  initial  conditions  for  the  equilibrium  and  frozen  flow  equations  are 
carried  out  until  sets  of  values  are  found  which  give  sonic  flow  at  the 
throat  of  the  same  nozzle  shape  as  was  determined  from  the  finite-rate 
solution. 

For  the  pressure  distribution  in  the  convergent  section,  we  selected 
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a  dependence  of  the  form  p  ■  pQ  e  ,  where  pQ  is  the  inlet  pressure 
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and  b  is  a  parameter  to  be  chosen.  The  inlet  conditions  are  not  stagnation 
conditions,  but  rather  correspond  to  a  large  but  finite  area  Aq  and  a 
small  but  finite  approach  Mach  number  Mq.  In  the  divergent  section,  the 
area  ratio  was  chosen  to  be  of  the  form 


A  -  A.  +  K(x  -  x.)' 


(5.1) 


where  A^  is  the  throat  area,  K  -  it  tan  9,  and  x*  the  distance  from  the 

r 

inlet  to  the  throat.  The  divergent  section  thus  asymptotes  to  a  conical  nozzle 

with  half >angle  9.  Figure  2  shows  the  nozzle  contour  which  resulted  from  this 

_2 

method  of  calculation  for  the  chosen  values  9  =  15°  and  b  *»  1  cm 

Then,  for  given  p(x),  the  governing  equations  for  non-equilibrium  flow 
can  be  written  as: 
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For  given  A(x),  Eqs.  (2.3),  (2.9)  and  (2.14)  remain  unchanged.  Equation 


(5.2)  is  replaced  by 
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and  instead  of  (5.2)  we  use 


d£  ,  _  £  dA  .  a  F  . 

dx  A  dx  u  2 


(5.5) 


The  equations  for  equilibrium  flow  are  obtained  by  setting  R  “  0, 


Tg  *  T,  and  replacing  (2.9)  by  (4.3). 
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For  fzozen  flow,  if  we  specify  p  ■  pQ  e  ,  the  equations  can  be 
integrated  and  one  obtains 
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and  hence,  from  the  condition  (dA/dx)  ■  0  at  x  ■  x*, 
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The  Mach  number  M  is  defined  as 
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If  H  is  small,  we  see  from  (5.9)  that  x.  depends  only  on  the 
parameter  b.  Also,  from  (5.10),  we  observe  that  (A*/Aq)  is  independent 
of  b  and  a  function  only  of  M^.  These  properties  serve  as  a  guide  in 
the  choices  of  Mq  and  b  for  a  desirable  inlet  shape  in  the  non-equilibrium 
flow  problem.  For  different  flow  conditions,  b  and  Mq  have  to  be  adjusted 
so  as  to  keep  the  convergent  portion  of  the  nozzle  the  same. 

In  the  actual  frozen  flow  calculations,  the  solutions  given  by  Eqs. 

(5.6)  -  (5.10)  were  not  used,  because  the  specification  P  “  PQ  e  was 

introduced  to  avoid  iteration  in  the  non-equilibrium  case.  The  frozen  flow 
solution  can  be  obtained,  of  course,  by  integrating  the  equations  for  a  given 
A(x)  determined  from  the  non-equilibrium  solution.  This  integration  can  be 
carried  out  in  closed  form,  with  the  result 


where 
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(5.13) 


(5.14) 


(5.15) 


and  T  is  given  by  the  equation  of  state.  It  was,  however,  found  more  con¬ 
venient  to  compute  the  frozen  flow  solution  by  numerical  integration  of  the 
differential  equations  than  by  numerical  evaluation  of  the  algebraic  relations 
(5.12)  -  (5.15). 

All  of  the  calculations  were  performed  on  an  IBM  7094  computer.  The 


Adams  four-point  integration  formula  was  used. 
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6.0  RESULTS  AND  DISCUSSION 

The  cases  which  were  computed  are  tabulated  in  Table  I.  The  inlet  con¬ 
ditions  in  all  cases  were  chosen  to  correspond  to  Saha  equilibrium,  although 
other  conditions  could  have  been  used.  If  the  gas  were  arc-heated  at  low 
pressure,  for  instance,  it  might  very  well  not  be  in  Saha  equilibrium.  The 
non-equilibrium  inlet  conditions  would  depend  in  a  complicated  way  on  the 
experimental  arrangement,  and  would  probably  have  to  be  measured  experiment¬ 
ally.  The  inlet  parameters  were  chosen  principally  to  illustrate  the  results 
of  the  analysis,  although  the  low  ionization  case  perhaps  is  somewhat  repre¬ 
sentative  of  arc-heated  wind  tunnel  performance.  Ideally,  it  would  have  been 
desirable  to  compute  enough  cases  to  permit  parametric  representation  of  the 
results  from  which  one  could  obtain  the  flow  properties  corresponding  to  any 
given  set  of  inlet  conditions  within  a  useful  range.  However,  the  numerical 
integrations  were  too  laborious  and  time-consuming  to  make  this  feasible,  and 
we  contented  ourselves  with  the  examination  of  essentially  two  cases,  repre¬ 
sentative  of  the  weakly  ionized  and  the  highly-ionized  regimes. 

The  argon  calculations  were  performed  primarily  to  illustrate  the 
influence  of  atomic  mass  on  the  recombination  coefficient,  and  it  was 
assumed  that  the  same  values  of  OL  and  S  as  used  for  hydrogen  would  be 
applicable  to  argon.  For  a,  this  assumption  should  be  quite  good  as  long 
as  the  electron  density  is  not  too  high.  Although  several  arguments  can  be 
advanced  to  support  this  point  of  view,  the  most  direct  ij  the  classical 
calculation  of  three-body  recombination  carried  out  by  Makin  and  Keck  [20]. 
These  results  are  in  agreement  with  the  quantum  mechanical  calculations  of 
Bates,  and  are  equally  valid  for  argon  as  for  hydrogen,  as  they  consider 
only  classical  electron  orbits  around  a  charged  core.  It  might  be  noted 
that  the  much  larger  recombination  rates  fi>r  argon  reported  by  Byron,  Bortz 
and  Russell  [19]  were  the  result  of  over-estimating  the  unknown  collisional 
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transition  probabilities,  and  it  is  now  generally  agreed  that  Bates’  hydrogen 
values  for  a  should  be  reasonably  accurate  for  argon  as  well  [22].  The 
values  of  S  for  argon  should  be  much  smaller  than  for  hydrogen,  because  of 
the  higher  first  excitation  potential,  but  S  is  only  important  in  the  region 
upstream  of  the  sonic  throat  when  conditions  are  near  Saha  equilibrium,  and 
thus  very  little  error  is  made  in  the  properties  in  the  supersonic  region, 
where  the  flow  becomes  highly  non-equilibrium.  The  radiated  power  R,  which 
enters  into  the  electron  energy  balance  given  by  Eq.  (2.10),  should  be  chosen 
so  that  the  net  energy  of  recombination  given  up  to  the  free  electrons  is  the 
same  for  argon  as  for  hydrogen  in  the  optically  thin  case.  We  thus  find  that 

d«e 

Rargon  *  ^hydrogen  dt 

In  the  actual  calculations  R  was  chosen  to  be  a  slightly  smaller  value, 

argon  ° 

so  that  the  fraction  of  recombination  energy  given  to  the  free  electrons  was 
about  16  percent  larger  than  that  given  above.  Thus  in  the  argon  calculations 
the  electron  temperature  is  a  little  higher  than  it  should  be,  but  the  differ¬ 
ence  is  quite  small,  and  less  than  the  uncertainties  of  the  model  of  recombina¬ 
tion  used. 

The  results  are  presented  in  Figs.  3  -  11.  The  curves  are  double-valued 
functions  of  A/A^  because  properties  in  both  the  subsonic  and  supersonic 
portions  of  the  nozzle  are  plotted.  The  gross  flow  properties,  p,  u,  p  and 
f  are  shown  in  Figs.  3-6.  We  observe  from  Fig.  3  that  for  f  *  0.62,  and 
9  *  15°  in  hydrogen,  the  density  variations  for  the  two  non-equilibrium  cases 
and  for  frozen  flow  are  essentially  indistinguishable.  The  density  falls  more 
rapidly  for  the  equilibrium  flow,  however.  Equivalent  information  is  presented 
in  Figs.  4a  and  4b,  where  the  flow  velocities  are  plotted.  Although  not  dis- 
ceraable  on  the  graphs,  the  velocity  for  a  *  0  is  slightly  greater  than 


that  for  a  =  <=.  The  highest  velocity  is  obtained  for  frozen  flow,  and  the 
lowest  for  equilibrium  flow. 

Computed  pressure  variations  in  the  nozzle  are  shown  in  Figs.  5a  and 
5b.  As  in  the  case  of  the  velocity  and  density,  not  much  difference  is  found 
between  the  non-equilibrium  and  frozen  flows.  The  pressure  in  equilibrium 
flow  remains  somewhat  higher  than  the  other  cases,  throughout  the  supersonic 
portion  of  the  nozzle,  for  f  «*  0.62.  For  the  low  initial  ionization, 
f  «  0.01,  the  pressure  variations  for  all  cases  are  practically  identical, 
as  might  have  been  expected. 

Figure  6  shows  the  variation  in  fractional  ionization  f,  for  initial 
condition  f  =  0.62.  The  effect  of  expansion  angle  is  what  one  expects; 
increasing  6  moves  the  flow  closer  to  the  frozen  limit,  and  vice  versa . 

The  non-equilibrium  finite-rate  flows  are  much  closer  to  frozen  flow  than 
to  equilibrium  flow.  The  actual  electron  number  densities  in  the  nozzle  are 
plotted  in  Figs.  7a,  7b  and  7c.  It  can  be  seen  that  for  the  low  initial 
ionization  case,  r  «  0.01,  the  equilibrium  flow  calculation  gives  an 
extremely  rapid  drop  in  .  One  notes  that  the  "sudden-freezing"  model  [2l] 
which  has  been  used  successfully  for  atomic  recombination  would  be  a  poor 
approximation  for  electron-ion  recombination,  as  has  been  previously  observed 
by  Bray  [15]  and  3yron,  et.al. [19]. 

Some  of  the  most  interesting  results  of  the  calculations  are  presented 
in  Figs.  8a,  8b,  8c  and  9,  where  the  electron  and  ion-atom  temperatures  are 
plotted.  In  Fig.  8a.  we  see  that  the  temperatures  for  the  non-equilibrium 
cases  in  hydrogen  are  closer  to  frozen  flow  than  to  equilibrium  flow,  and 
the  effect  of  increasing  the  expansion  angle  is  to  move  the  flow  closer  to 
the  frozen  limit,  as  one  might  have  anticipated.  The  temperature  level  for 
o  **  oo  is  higher  than  for  C  ■  0,  as  is  the  temperature  difference  Tg 
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We  observed  that  more  recombination  occurs  for  a  -  0,  so  that  more  recom¬ 
bination  energy  is  liberated  for  this  case.  However,  this  recombination 
energy  ic  largely  lost  by  radiation,  with  the  result  that  the  gas  temperature 
is  lower  for  the  optically  thin  case  than  for  the  optically  thick  case.  The 
temperature  difference  Tg  -  T  is  larger  for  the  optically  thick  case  because 

the  effectiveness  of  energy  transfer  between  electrons  and  ions  is  proportional 
-3/2 

to  T  '  ,  and  Tg  is  larger  for  the  optically  thick  case.  Thus  radiation 
power  loss  serves  both  to  increase  the  amount  of  recombination  and  cool  tbs 
gas,  while  decreasing  the  temperature  difference  T£  -  T.  Similar  comments 
apply  to  the  data  plotted  in  Fig.  8b.  Although  the  gross  flow  variables  are 
practically  unaffected  by  recombination,  for  f  «=*  0.01  we  observe  that  there 
is  a  significant  effect  on  Te  due  to  radiation  absorption. 

The  temperature  behavior  of  argon  is  shown  in  Fig.  8c.  Here  we  observe 
that  for  fQ  =  0.01,  the  gas  temperature  T  is  hardly  affected  by  the 
presence  or  absence  of  recombination.  However,  Tg  at  the  exit  of  the 
nozzle  (A/A+  «  100)  is  almost  an  order  of  magnitude  greater  than  T.  As 
in  the  case  of  hydrogen,  the  difference  (Tfi  -  T)  is  greater  for  the  optically 
thick  case  than  for  the  optically  thin  case.  For  simila;  onditions,  this 
temperature  difference  (Tg  -  T)  is  larger  in  argon  than  in  hydrogen.  This 
is  because  the  energy  exchange  rate  between  electrons  and  massive  particles 
is  inversely  proportional  to  the  atomic  mass  m  ,  as  seen  from  Eqs.  (2.11) 
and  (2.12).  The  data  for  (Tg  -  T)  have  been  replotted  in  Fig.  9  to  show 
this  effect  more  clearly. 

Another  comparison  between  the  behavior  of  argon  and  hydrogen,  in  terms 
of  the  fractional  amount  of  recombination,  is  shown  in  Figs.  10a  and  10b.  Two 
competing  effects  are  involved  in  determining  the  variation  in  A  ■  (fQ  -  f)/f  « 
Since  the  flow  velocities  are  lower  in  argon  than  in  hydrogea,  the  residence 
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time  in  the  nozzle  is  longer  for  the  argon,  and  this  tends  to  increase  the 

net  recombination.  However,  for  reasons  just  discussed,  the  temperature 

difference  (Tg  -  T)  is  greater  in  argon,  and  the  recombination  rate  is  a 

strong  inverse  function  of  T  .  For  the  case  of  high  initial  ionization, 

e 

fQ  *  0.62,  the  residence -time  effect  predominates,  and  A^  >  A^  at  large 

A/A*,  as  seen  in  Fig.  10a.  The  converse  occurs  for  low  initial  ionization, 

f  *  0.01,  as  seen  in  Fig.  10b. 
o 

The  Mach  numbers  obtained  for  the  various  cases  are  plotted  in  Figs. 
11a  and  lib.  In  general,  the  behavior  is  what  might  have  been  expected. 


7.0  EFFECT  OF  HEAT  CONDUCTION 

It  has  been  recognized  [14,15]  that  the  presence  of  high  temperature 
electrons,  even  in  small  mole-fractional  amounts,  can  increase  the  thermal 
conductivity  of  a  gas  markedly.  In  a  rapid  expansion  through  a  nozzle,  large 
axial  temperature  gradients  are  produced,  and  it.  is  appropriate  therefore  to 
examine  the  effect  of  axial  heat  conduction  on  the  energy  balance. 

The  importance  of  axial  heat  conduction  can  be  assessed  by  examining 
the  value  of  the  quantity 


which  is  the  ratio  of  the  net  electron  enthalpy  flux  and  the  electron  hsat 
flux  out  of  a  control  volume.  Conduction  effects  will  be  negligible  if 
C  »  1. 

For  the  electron  thermal  conductivity  kg,  we  use  the  value  given  by 
Spitzer  I 143  (pp.  144-145) 


k  =»  €  5_  kT 
e  T  L 
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where 
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In  our  case,  e  «  0.419  and  5^  =  0.225.  For  a  Lorentz  gas,  f  «  1 
5T=  !. 
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Similar  numerical  values  are  obtained  in  argon.  Therefore,  in  all  our 
examples  Involving  finite  rates  and  frozen  flow,  the  neglect  of  heat  con¬ 
duction  seems  justified.  For  equilibrium  flow,  however,  it  appears  that 
the  conduction  term  must  be  included.  One  may  conclude  that  in  general  the 
heat  conduction  term  may  be  neglected  except  at  very  Mgh  electron  temper¬ 
atures  or  for  very  small  degree  of  ionization.  In  the  latter  case,  conduction 
effects  would  act  to  keep  the  electron  temperature  nearer  to  its  initial 
value  ahead  of  the  expansion.  These  conclusions  are  essentially  the  same  as 
those  reached  by  Bray  [153. 


O  1 
JL  • 


8.0  CONCLUDING  REMARKS 

We  have  not  attempted  any  detailed  comparisons  between  our  argon  results 
and  those  of  Bray  and  Wilson  [8,15],  because  the  initial  conditions  taken 

by  these  authors  differ  from  ours.  However,  one  may  remark  that  for  the  most 
part,  their  results  and  ours  are  qualitatively  similar.  The  principal  feature 
of  our  results  which  we  believe  is  new  is  a  more  complete  and  accurate  evalu¬ 
ation  of  the  behavior  of  the  electron  and  gas  temperatures,  and  of  the  effects 
of  radiation  loss  (as  predicted  for  the  optically  thin  case,  a  =  0)  on  the 
overall  recombination  process.  These  effects  may  be  summarized  as  follows: 

In  comparing  the  two  cases  a  **  0  and  o  *  at  a  common  station  A/A^. 
in  the  nozzle,  one  finds  that  radiation  loss  (a)  produces  more  cooling  of 
the  gas,  (b)  increases  the  net  recombination,  (c)  decreases  the  temperature 
difference  (Tg  -  T)  and  (d)  increases  the  expansion  pressure  ratio. 

The  effect  of  decreasing  the  expansion  angle  of  the  nozzle  is  to  increase 
T  and  T,  and  increase  the  net  recombination  A,  but  to  decrease  the 
temperature  difference  (T£  -  T).  All  of  these  effects  are  traceable  to 
the  increase  in  residence  time  of  the  fluid  in  the  nozzle. 

For  the  case  considered  with  initial  ionization  fraction  f  *  0.62, 

o 

recombination  effects  alter  the  gross  flow  variables  in  the  nozzle  measurably 
though  not  dramatically.  For  the  low  initial  ionization  case,  fQ  “  0.01, 
the  gross  flow  variables  are  essentially  unaltered  by  the  presence  or  absence 
of  recombination,  and  the  only  significant  result  of  the  recombination  is  to 
maintain  (Tfi  -  T)  >  0. 

On  comparing  the  argon  and  hydrogen  flows,  one  finds  that  the  effect  of 
the  larger  atomic  mass  of  the  argon  is  to  diminish  the  effectiveness  of 
energy  exchange  between  the  electrons  and  massive  particles,  resulting  in  a 
larger  t.  aperature  difference  (Tg  -  T)  in  argon.  For  the  case  fQ  **  0.62, 


the  fractional  amount  of  recombination  A  downstream  in  the  supersonic  por¬ 
tion  of  the  nozzle  was  greater  in  argon  than  in  hydrogen,  because  of  the 
longer  residence  time  of  the  argon  in  the  nozzle.  But  for  f  *  0.01,  A 
was  greater  in  hydrogen,  because  the  higher  argon  electron  temperature 
reduces  the  recombination  rate. 

Electron  thermal  conductivity  effects  were  found  to  be  unimportant 
for  the  frozen  and  finite-recombination-rate  flows,  but  should  be  included 
in  equilibrium  flow  calculations,  at  least  in  the  supersonic  port'on  of 


the  nozzle. 


23. 


APT  3 IX  A 

The  rate  constants  and  radiated  power  factors  computed  by  Bates,  et.al. 

are  reported  in  tabular  form  at  rather  widely  spaced  intervals  in  and 

Te>  For  the  numerical  calculations  of  this  report  it  was  convenient  to 

obtain  an  approximate  analytical  representation  of  these  results.  Since 

the  Bates,  et^.al.  results  have  a  rather  broad  range  of  useful  application, 

it  seems  worthwhile  to  report  our  analytical  approximations. 

It  was  found  that  the  recombination  coefficient  a  for  a  *  0  could 
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be  represented  adequately  over  range  500°K  <  8000°K,10  <  Ne  <  10  by 

the  double  series 


ioSio  « 


I 


(A.  1) 


The  twelve  coefficients  A^  were  found  by  using  Bates'  tabular  values  and 
solving  the  resulting  matrix  equation;  we  obtained 


Aoo 

-  13.825 

A01  " 

5.212 

A02  3 

-  0.7959 

A  ** 

10 

0.05299 

A11  " 

0.7049 

A12  " 

0.06981 

> 

(S3 

0 

a 

-  0.14655 

A21  ° 

0.6805 

A22 

-  0.458 

A30 

0.03 

A31  - 

-  0.1058 

A32 

0.07253 

Similarly,  for  the  ionization  rate  S  for  0  «  0  we  write 

v’v1  r  /  »  \  i1  /  T  \j/2 

logios  "  L  L  Ll0UI^)J  (4000)  (A 

i-0  j-0 

12  -3 

In  the  range  Tg  >  4000°K,  Ng  >10  cm  ,  we  find  the  nine  coefficients 


to  be 


Boo 

S  - 

62.6086 

BQ1  -  51.004 

*02  ' 

-  12.8399 

B10 

m 

2.1364 

Bn  -  -  1.893 

012 

0.4810 

B20 

-  - 

0.2685 

B21  -  0.4842 

B22  " 

-  0.1542 

For 

the  power  radiated,  R 

,  we  write, 

for  ff*«0 

R  - 

R  +  N  1 

R, 

(A.  3) 

o  e 

1 

where 

V 

g 

corresponds  to  the  quantity 

tabulated  by  Bates,  jat.aT ,  and 

R1 

corresponds  to  their  °R^. 

We  write 

3 

2 

r  N 

,l08l»  1^2  . 

i 

j/2 

lo810  “o 

'  I 

Iv. 

(t; 

\ 

) 

i-0 

j-0 

(  (A. 4) 

3 

V"1 

2 

Ne  1 

1  ,  Te 

V,  4000 

j/2 

\ 

l°810  R! 

‘  I 

.l08l°  IJw  j 

>  ; 

i*0 

j«0 

and 

obtain 

coo 

■  - 

0.5088 

C01 

5.2306 

C02  *  • 

•  0.8142 

C10 

m 

2.1192 

C11 

0.3779 

C!2  ’ 

0.3891 

C20 

0.2545 

C21 

1.0944 

C22  ‘  ' 

■  0.8628 

C30 

m 

0.0559 

C31 

-  0.2108 

C32 

0.1756 

Doo 

48.924 

D01 

38.89 

D02  ’  ' 

10.07 

D10 

*  - 

2.932 

D11 

6.681 

D12  “  ' 

2.753 

D20 

m 

3.352 

°21 

-  5.684 

°22  ' 

2.337 

°30 

0.6826 

D31 

1.158 

8  32  ‘  ' 

0.4774 

The  range  of  applicability  of  R  is  approximately  500°  <  Tg  <  16,000°, 
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10A  <  N  <  10  . 

e 


co  « 


For  the  optically  thick  case,  9  ■  »,  it  is  found  that  the  dependence 


of  a  on  Ng  can  be  written  explicitly  as 


a  =  A  +  K  N 
«  R  t  e 


On  fitting  Bates1  tabulated  values,  we  find 


(A.  5) 


lO*10  ar 


182.333  -  1482.19  (  1|22  )  '  +  4455.599  (  1222  ^  ^ 

37^.5845  (  1220 

,  s/a 

-  2530.238  ^  ^  +  450.16  (  .  (A. 6) 

165.9957  -  1440.; 


7  1000  \3/4 
-  687u.2o  i  — —  ) 


log10  Kt 


.  (  1000  \ 

+  450.18  ( 1222  f2  . 


-  6408 


7328  (  122£  y  +  4245.2524  (  1225  j  ^ 

‘•“KT/' 


+  5299.7876  f  ) 
<=  -  '  e  ' 

-  2285.3684  ^  ^2,  ^  +  402.2243  ^ 

»n  is  valid  over  the  ranee  250°  <  T  < 


1222  f1  .  (a.7) 

3  ' 


This  representation  is  valid  over  the  range  250°  <  T  <  16,000°.  The  ioniza¬ 
tion  rate  for  9  »  <*>  can  be  expressed  by  the  simple  relation 


10 


-23.94 


^  1000  ) 


11.03 


(A.  8) 


in  the  range  8000°  <  T  <  32,000°. 


26. 


If  we  use  this  constant  to  calculate  the  quantity  given  by  Eq.  (2.1),  we 

find  that  it  agrees  well  with  the  value  given  in  Ref.  7  for  -  1  eV. 

The  collision  probability  P£  for  argon  has  been  measured  by  Brode  I 18], 
who  found  P£  -  3  for  kTg  ■  1  eV.  Using  this  value,  we  obtain  for  argon 

_l/2  ,  v  in-23  1/2  2 

B  '  ■  7  X  10  erg  '  cm 


(B.6) 
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cipal  effect  of  the  recombination  is  to  produce  an  elevated  electron  temperature. 

On  comparing  the  hydrogen  and  argon  flows,  one  finds  that  the  effect  of  the 
larger  atomic  mass  of  the  argon  is  to  diminish  the  effectiveness  of  electron-atom 
energy  exchange,  resulting  in  a  larger  difference  between  electron  and  atom  tempera  :ures 
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